The Distribution of Dark Matter in Galaxies: 
the Core Radius Issue 



Paolo Salucci 1 

SISSA, Via Beirut, 2-4 1-34014 Trieste, Ita;y, salucci@sissa.it 



Abstract. I review the up-to-date status on the properties of the Dark Matter density 
distribution around Galaxies. The rotation curves of spirals all conform to a same 
Universal profile which can be uniquely decomposed as the sum of an exponential thin 
stellar disk and a dark halo with a flat density core. From dwarfs to giants galaxies, the 
halos embedding the stellar component feature a constant density region of size ro and 
value po, which are inversely correlated. The fine structure of dark halos in the region 
of the stellar disk has been derived for a number of low-luminosity disk galaxies: the 
halo circular velocity increases almost linearly with radius out to the edge of the stellar 
disk, implying, up there, an almost constant dark matter density. This sets a serious 
discrepancy between the cuspy density distribution predicted by N-body simulations 
of ylCDM cosmology, and those actually detected around galaxies. 

The small scatter around the Fundamental Plane (FP) of elliptical galaxies con- 
straints the distribution of dark and luminous matter in these systems. The measured 
central velocity dispersion <tq in the FP is linked to both photometric and dynamical 
properties of luminous and dark matter. As a consequence, the well-known features of 
the FP imply that, inside the effective radius R s , the stellar spheroid must dominate 
over the dark matter, in contrast with ylCDM predictions. 



1 Introduction 

Rotation curves (RC's) of disk galaxies are the best probe for dark matter (DM) 
on galactic scale. Notwithstanding the impressive amount of knowledge gathered 
in the past 20 years, only very recently we start to shed light on crucial aspects 
of the mass distribution of Dark Matter in galaxies, including its radial density 
profile, and its Universality. 

On a cosmological side, high-resolution N-body simulations have shown that 
cold dark matter (CDM) halos achieve a specific equilibrium density profile [13 
hereafter NFW, 5, 8, 12, 9] characterized by one free parameter, e.g. the halo 
mass. In the innermost region, the DM halo density shows an average profile 
which is characterized by a power-law cusp p ~ r -7 , with 7=1 — 1.5 [13, 12,2]. 
In detail, CDM halos have: 

PNFWO) = PS t— (1) 

(r/r s )(l + r/r s y 

where r s and p s are respectively the characteristic inner radius and density. Let 
us define r v ; r as the radius within which the mean density is A v 
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Fig. 1. Synthetic rotation curves (filled circles with error bars) and the Universal Ro- 
tation Curve (solid line). The separate dark/luminous contributions are indicated by 
a dotted line (disk) and a dashed line (halo). 



universal density p m at the halo formation redshift, the associated virial mass as 
M v ; r and the halo velocity as V v - U = GM v ; r /r v i r . In the "concordance" /1CDM 
scenario: fi m — 0.3, Qa = 0.7 and h = 0.7, so that Z\ v ; r ~ 340 at z ~ 0. Let 
us set c = r v i r /r s , and from simulations c ~ 21(M t , ir /10 11 ) -013 , then in CDM 
framework circular velocity of dark halos VnfwM depends only on their virial 
masses and takes the form: 

^NFw( r ) = ^ir~^FT ~ " ( 2 ) 

where x = r/r s and A{x) = ln(l + x) — x/(l + x). 

From the observational point of view, only recently the difficulties in deriving 
the internal structure of halos from available kinematics have been overcome. 
This has been done i) by means of a careful study of the Universal Rotation 
Curve [16] built out of 1000 individual RC's, ii) by adopting an halo velocity 
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Fig. 2. a vs. (5 and f3 vs. V op t. 



profile that, out to r opt , is neutral with respect to various different galaxy mass 
models: 

KuBci*) = Vo P t (1 -/?)(! + « 2 ) (a / + 2 fl2) (3) 

where x = r/r opt , a the halo core radius in units of r op t and (3 = (Vd.uRc) /Vopt) 2 
at Ropt- Then, by varying (3 and a, Vft.,[/,Rc can reproduce the maximum-disk, 
the solid-body, the no-halo, the all-halo, the NFW mass models, (e.g. CDM 
halos with concentration parameter c — 5 and r s = r op t are well fit by (3) with 
a ~ 0.33) Hi) by means of a number of high-quality high-resolution individual 
RC's [1] leading to trustworthy mass distributions. Let us define hereafter Td as 
the disk scale-length and V opt = V(r opt ). 

2 Dark Halos Properties from the Universal Rotation 
Curve 

Let us remind the observational framework: a) the mass in spirals is distributed 
according to the Inner Baryon Dominance (IBD) regime [16]: there exists a char- 
acteristic transition radius vjbd — 2rd(V op t/220 km/s) 1 ' 2 for which, at r < vibd, 
the luminous matter completely accounts for the gravitating mass, whereas, at 
r > riBD, the dark matter shows up in the kinematics and rapidly becomes the 
dominant mass component [20, 18, 1]. Then, although dark halos extend down 
to the galaxy centers, it is only for r > rjg£> that they give non-negligible con- 
tributions to the circular velocity, b) DM is distributed in very differently way 
with respect to the baryons [16, 6], and c) the HI contribution to the circular 
velocity, at r < r opU is small [e.g. 17]. 

Persic, Salucci and Stel [16] have derived from ~ 20000 velocity measure- 
ments, relative to ~ 900 rotation curves, V syn ( y—^ ; Mi ) , the synthetic rotation 
velocities of spirals binned in intervals of magnitudes. Each individual RC's (see 
Fig. (1)) shows a variance, with respect to synthetic curves of the corresponding 
magnitude, smaller than observational errors: spirals sweep a very narrow locus 
in the RC- profile/amplitude/ luminosity space. Thus, as regard the average 
main properties of the DM distribution, eq (3) is equivalent to a large sample of 
individual objects. 
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Fig. 3. up) Central halo density po (in g/cm 3 ) vs. disk mass (in solar units) for normal 
spirals (filled circles); bottom) central density vs. core radii (in kpc) for normal spirals 
(see [19]). 

The whole set of synthetic RC's define the Universal Rotation Curve (URC) 
that we represent analytically with the sum of two terms: a) the standard expo- 
nential thin disk term: 

VluRcix) = 1.28 V a 2 pt x 2 (I K - hK^x (4) 

and the spherical halo term given by (3). The data (i.e. the synthetic curves 
V„yn) select the actual mass model: by setting V§ RC (x) = V 2 URC (x, f3,a) + 

urc( x i 0) with a and /3 as free parameters, an extremely good fit occurs 
when: (3 = /3(log V opt ) and a — a((3) as plotted in Fig. (2): the URC reproduces 
V syn (r) up to its rms (i.e. within 2%). Moreover, at a fixed luminosity, the la 
fitting uncertainties for a and (3 are lesser than 20%. 

The emerging scenario is the following: inside r opt smaller objects have larger 
dark-to-stellar matter : M*/M vir ~ 0.2 (M*/2 x 10 n M l3 ) j (j - 0.75) [20]) 
and within each galaxy the dark mass increases with radius with a power-law 
exponent between 2 and 3. 

This evidence calls for a cored Dark Matter density [3, 4, 1]. Then, we are 
allowed to pass from the " neutral" distribution of DM of eq (3) in which a " core" 
may appear in the velocity profile (i.e. in a quantity which is directly measured), 
to the much more specific mass distribution given by the Burkert density profile 
that forces a core radius into the NFW profile. 
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with p and r free parameters (the central DM density and the core radius). 
Then M B (r) = 4 M {ln(l + r/r ) - arctan(r/r ) + 0.51n[l + (r/r ) 2 ]} with 
Mo ~ 1.6 po Tq the dark mass within rrj. The halo contribution to the circular 
velocity is: V^ B (r) = G M B (r)/r. 

The Disk + Burkcrt halo model applied to the synthetic rotation curves leads 
to core parameters tq, pq strong correlated and linked to the disk mass : dark 
halos behave as an 1-parameter family, completely specified by e.g. their central 
density po (see Fig. (3) ). 

These relationships imply that the densest halos harbor the least massive 
disks (see Fig. 3) while the feature of the curvature at the highest masses/lowest 
densities in the po vs. ro relationship may be related to the existence of an upper 
mass limit in M v i r at few x 10 12 M Q . 

2.1 Testing CDM with the URC 

From the analysis of the URC one concludes that dark halos arc not kincmatically 
cold structures, but "warm" regions of sizes r oc p^ 1,5 with r <~ 4 — 7 r^. Then, 
the boundary of the core region is beyond the region of the stellar disk. There 
is no evidence that the DM density converges to a r 1 (or a steeper) regime, as 
dictated by CDM scenario. 

3 Dark Matter Properties from Individual Rotation 
Curves 

Although deriving halo densities from individual RC's is certainly complicated, 
the belief according to which one is bound to get ambiguous halo mass modeling 
[as claimed in some work] is not always correct. In fact, this is true only for ro- 
tation curves of low spatial resolution, i.e. those with less than ~ 3 measures per 
exponential disk length-scale, as most of HI rotation curves. In this case, since 
the parameters of the galaxy structure are crucially sensitive to the shape of the 
rotation curve in the region < r < rj, there are no sufficient data to constrain 
the mass model. In the case of high-quality optical RC's tens of independent 
measurements solve the problem. Moreover, since the dark component can be 
better traced when the disk contributes to the dynamics in a modest way, a con- 
venient strategy leads to investigate DM-dominated objects, like dwarf and low 
surface brightness (LSB) galaxies. For the latter [e.g. 7, 11, 3, 4, 9, 10, 21, 23] the 
results are not definitive in that they unfortunately are intrinsically uncertain, 
due to the limited quality of their kinematics. 

Since the observed RC's have Universal features and most of the properties of 
cosmological halos are also claimed universal, an useful strategy is to investigate a 
restricted number of high-quality optical rotation curves of low luminosity late- 
type spirals, with /-band absolute magnitudes —21.4 < Mj < —20.0 and 100 < 
V op t < 170 km s _1 . Objects in this luminosity/velocity range are DM dominated 
[e.g. 20] but their optical RC's, measured at a typical angular resolution of 2", 
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Fig. 4. B+ Disk models (thick solid line) (points with errorbars). Thin solid lines 
represent the disk and halo contributions. The maximum and minimum disk solutions 
(dashed lines) provide the uncertainties see [1] 



have the excellent spatial resolution of ~ 100(D/10 Mpc) pc and ndata ~ r opt /w 
independent measurements. For nearby galaxies: w « and ndata > 25. 

In [1] we extracted, from the 'excellent' subsample of 80 rotation curves of 
[15], the best 9 rotation curves. These RC's (as any rotation curve that candi- 
dates itself to yield something crucial on the DM distribution) trace properly 
the gravitational potential since satisfy the following quality requirements: 1 ) 
data extend at least out to the optical radius, 2) they are smooth and symmet- 
ric, 3) they have small internal rms, 4) they have high spatial resolution and a 
homogeneous radial data coverage of 30 — 100 data points between the two arms. 
Each RC has 7 — 15 velocity points inside r opt , each one being the average of 
2 — 6 independent data. The RC's spatial resolution is better than 1/20 r opt , the 
velocity rms is about 3% and the RC's logarithmic derivative is generally known 
within about 0.05. 

We model the mass distribution as the sum of two components: a stellar 
disk and a spherical dark halo, therefore: V 2 (observed) = V 2 (disk) + V^(halo). 
Light traces the stellar mass via the radially constant mass-to-light ratio. We 
neglect the gas contribution V gas (r) since in normal spirals it is small [17, Fig. 
4.13]: I3 gas = (V 2 as /V 2 ) ropt ~ 0.1. Incidentally, this is not the case for dwarfs 



The Distribution of Dark Matter in Galaxies: the Core Radius Issue 



7 



and LSBs: most of their kinematics is affected in different ways by the HI disk 
gravitational pull. The disk contribution to circular velocity is given by (4), while 
the dark halo contribution by (3). Finally we normalize (at r opt ) the velocity 
model {yj + V, 2 ) 1 / 2 to the observed rotation speed V opt . 

For each galaxy, we determine the values of the parameters (3 and a by means 
of a x 2 -minimization fit to the observed rotation curves: 

VLdei (r;/3,a) = V 2 d (r ; (3) + V* (r;/3,a) (6) 

A central role in performing the mass decomposition is played by the derivative 
of the velocity field dV/dr. It has been shown [e.g. 14] that by taking into account 
the logarithmic gradient of the circular velocity field defined as: V(r) = d ^^7-^ 
one best retrieves the crucial information stored in the shape of the rotation 
curve. Then, we set the % 2 statistics as the sum of x 2 's evaluated on velocities 
and on logarithmic gradients. In detail, by setting Xv = S"=i Vi ~ Vmo ^y^ Ti ' l3,a " > 

and xv = E7=i V(r ' ) ' V ^r (r ' :/? ' a) - we minimize the quantity 

Xtot = Xv + Xv ( 7 ) 

to get the mass model. 

Let me point out that any claim of " different mass models" that all would 
account for a certain rotation curve, does instead originate from the (low) quality 
of the latter that does not allow a reasonably accurate derivation of V(r). 

The best-fit models parameters for the "neutral" distribution of (3) are 
shown in Fig. 4. The disk-contribution (3 and the halo core radius a span a 
range from 0.1 to 0.5 and from 0.8 to 2.5, respectively. They are pretty well 
constrained in a small and continuous region of the (a, (3) space. It is obvious 
that the halo curves are increasing almost linearly, out to the last data point. 
Remarkably, we find that the size of the halo density core is always greater than 
the disk characteristic scale-length and it can extend beyond the disk edge 
(and the region investigated). 



3.1 Testing CDM 

Let us assume for the dark halos the NFW functional form given by (1) and fit 
the 9 RC's leaving c and r s as free independent parameters, although N-body 
simulations and semi-analytic investigations indicate that they correlate in order 
to increase the chance of a good fit. It must be assumed, however, for these 9 
test objects, a conservative halo mass upper limit of 2 x 10 12 M Q . The fits to the 
data are shown in Figs. (5) and (6), together with the URC fits: for seven out of 
nine objects the NFW models are unacceptably worse than the URC solutions. 
Moreover, the resulting CDM disk /-band mass-to-light ratios turn out to be in 
some cases ~ 0.01 solar units, i.e. unacceptably low in the /-band. 

We definitely conclude that there is no shortage of dark halos around objects 
with a trustable rotation curve, that show a density distribution inconsistent 
with that predicted by collision-less CDM. 
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Fig. 5. NFW best-fits solid lines of the rotation curves (filled circles) compared with 
the URC halo + Disk fits (dashed lines). The x 2 values are also indicated. 




Fig. 6. left) The density of the dark halo of 116-G12, right) The CDM prediction 



4 Halos around Ellipticals 

The very low scatter that ellipticals show lying on the Fundamental Plane is a 
kincmatical feature can be used to investigate their Dark Matter distribution. 
The central velocity dispersion cfq is the key quantity, in that it strongly de- 




Fig. 7. The CDM surface in the log-space (R e , L r , ao) compared to sample data 
(points). R e in kpc, L r in L t q and ao in km/s. 



pends on the mass distribution of both stellar and dark matter [34]. A sample 
of 221 E/SO galaxies with good photometric and spectroscopic data has been 
used to investigate this point. This sample defines the standard "empirical" FP 
in the coordinate space of {loga^, logR e , logL r ). The galaxies have a average 
distance of 0.084 dex with respect to the plane, of the order of the measurement 
uncertainties, estimated as large as ±0.05 dex. 

It is easy to show, for the NFW distribution of ylCDM cosmology, that CTo 
is a specific function of the total virial mass and the stellar mass-to-light ratio, 
with no other free parameters. For the Burkert (1995) density distribution, <7o 
depends instead on the stellar mass-to-light ratio, the value of the core radius 
ro and the fraction of dark mass inside R e . 

It is shown that the existence of a tight FP relating the above quantities 
implies that ellipticals are largely dominated, within R ei by the stellar spheroid, 
independently of the actual DM distribution. However, ylCDM predicts large 
amounts of dark matter inside R e in view of the cuspy density distribution of 
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CDM halos: as results in this framework the ellipticals would lie on a it curved 
surface, inconsistent with the observed plane (see Fig. (7)). 



l0gL r 




Fig. 8. The plane from Spheroid+Burkert halo model. 

The Burkert density distribution, in which substantial amounts dark matter 
can be placed outside R e , leads to a relationship that is a plane resulting in 
perfect agreement with the observed one (see Fig (8)). This implies a dark- 
to-luminous mass fraction within the effective radius of about 0.3 ± 0.2 and a 
luminosity dependence of the spheroidal mass-to-light ratio: M sp h/L r — (5.3 ± 
0.1)(L r /L» r )°- 21±0 - 03 , in Gunn-r band. Moreover, as a firm constraint, we can 
state tq > R e . 



5 More Support for Core Radii. 

The Trieste group has provided a crucial evidence on the "core radii" issue, 
however, results from other investigations are also very important and must be 
considered, in that, alongside with those referred in the previous sections, build 
a formidable case for the existence of constant density cores at the centers of the 
dark halos surrounding spiral galaxies. 
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Table 1. Extra Evidence 



Author(s) 
Ortwin et al 
Fuchs 

Swaters et al . 
Weldrake et al. 
Bottema & 
Bolatto et al. 
de Blok et al. 
de Blok et al. 
Stil & Israel 
Dutton et al. 
Marchesini et al. 
Fraternali et al. 



Paper 

AJ 121, 1936 
astro-ph/02 12485 
ApJ 583, 732 
MN 340, 12 

Verheijen A&A 388, 793 
ApJ 565, 238 
AJ 122, 2396 
ApJL 552, 23 
A&A 392, 473 
astro- ph/0310001 
ApJ. 575, 801 
AJ 123, 3124 



6 The Intriguing Evidence from Dark Matter Halos 

From observations the dark halos around galaxies emerge as an one-parameter 
family; the order parameter (either the central density or the core radius) corre- 
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Fig. 10. DDO 47 rotation curve vs different mass models. The continous line is the 
sum of the stellar (red), the gas (yellow) and the Burkert halo distributions (magenta). 
The dashed line represents is the NFW distribution. 



lates with the stellar mass. However, the global structural properties of the dark 
halo, like the virial radius or the virial mass could be extrapolated only empir- 
ically, because we do not know how to prolongate, outside the region mapped 
by data and in a theoretical way, the DM profile. In fact, inside the stellar 
disk region, the halo RCs are determined by physical parameters, the central 
core density and the core radius, that have no counterpart in the gravitational 
instability/hierarchical clustering picture. 

Important relationships among physical quantities still emerge from the mass 
modelling: in Fig. 9 we show the dark-to-luminous mass ratio as function of the 
normalized radius and of the total disk mass. The surface has been obtained by 
adopting the correlations between the halo and the disk parameters derived in 
[19]. Therefore the dark-to-stellar mass ratio, at fixed fraction of disk length- 
scales, increases as the total disk mass decreases; for example: in the range 
1 < r/rd < 3 it raises by 20% for massive disks (Md = 10 12 M©) while it raises 
by 220% for smaller disks (M d = 1O 9 M ). 

Two conclusive statements can be drawn: dark matter halos have an inner 
constant-density region, whose size exceeds the stellar disk length-scale. Second, 
there is no evidence that dark halos converge, at large radii, to a p ~ r~ 2 (or 
steeper) profile. 

The existence of a region of "constant" density as wide as the stellar disk 
is hardly explained within current theories of galaxy formation. A number of 
different solutions have been proposed to solve this problem [e.g. 30, 31, 32, 33]. 
Before commenting on them let us stress that any solution of the " core radius" 
issue must account for all the intriguing halo properties described in this review. 
Let us point out that we review the several tenths of attempts for a solution 



The Distribution of Dark Matter in Galaxies: the Core Radius Issue 13 

to the "core radii" issue will be reviewed elsewhere. He we just classify them in 
families: 

1- Dark Matter "interacted" with itself or with baryons. Original "cuspy" 
halos have been smoothed in this way. 

2- Dark Matter has a different power spectrum/perturbations evolution than 
the current Standard Picture and this has "produced" the cores. 

3- Dark Matter is a " field" , that naturally mimicks the effects of a cored halo 
of particles. 

4- Within the standard yICDM scenario, N-Body simulations have failed, for 
a variety of reasons, to discover this intrinsic and real feature. 

7 Prologue 

Let me finish where I started: i.e. by showing a test case for the existence of a 
core in the density distribution of the dark halo around a galaxy. DDO 47, see 
Fig (10) has a rotation curve that increases linearly from the first data point, at 
300 pc, up to the last one, at 5 kpc well beyond the stellar disk edge. The RC 
profile implies the presence of a dominating (dark) halo with an (approximately) 
constant density out to the last measured point, prior any mass modelling. 
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